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ABSTRACT 


This publication presents a novel technique for simultaneously detecting data and estimating the pa- 
rameters of a received carrier signal phase modulated by unknown data and experiencing very high Doppler, 
Doppler rate, etc. Such a situation arises, for example, in the case of Global Positioning Systems (GPS) 
where the signal parameters are directly related to the position, velocity and acceleration of the GPS receiver. 

The proposed scheme is based upon first estimating the received signal local (data dependent) parameters 
over two consecutive bit periods, followed by the detection of a possible jump in these parameters. The 
presence of the detected jump signifies a data transition which is then removed from the received signal. 
This effectively demodulated signal is then processed to provide the estimates of global (data independent) 
parameters of the signal related to the position, velocity, etc. of the receiver. 

One of the key features of the proposed algorithm is the introduction of two different schemes which 
can provide an improvement of up to 3 dB over the conventional implementation of Kalman filter as applied 
to phase and frequency estimation, under low to medium signal-to-noise ratio conditions. One scheme is 
based upon reprocessing (cycling) the measurements over an optimally selected interval while the alternative 
scheme proposes an adaptive Hilbert transform technique. The overall complexity of the proposed algorithm 
is about three times the complexity of a single third order Kalman filter. 
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1. INTRODUCTION 


The problem of estimating the parameters of a received quasi-sinusoidal signal in the presence of noise 
has received considerable attention in the literature [1-12]. However, for the case when the received carrier 
is modulated by unknown data and simultaneously experiences considerably high Doppler and Doppler rate, 
the research reported in the published literature is somewhat limited. This situation arises, for example, in 
the case of Global Positioning Systems (GPS) receivers exhibiting very high dynamics [7, 12]. 

In a previous publication [7], an estimator structure based on the maximum likelihood estimation (MLE) 
of code delay and Doppler frequency over a single data bit period has been proposed and analyzed for the 
GPS application. The “pseudo” estimates over different bit intervals are combined by a Kalman filter to 
provide tracking of Doppler frequency. By limiting the primary (MLE) estimation period to less than one 
data bit period, the problem of detecting the data bits is bypassed. However, perhaps due to such a limitation 
and also due to the high frequency rate involved (not explicitly estimated by the MLE), a threshold of about 
30 dB-Hz in terms of the received carrier power to noise power spectral density ratio (P/No) was obtained. 
Due to the lack of knowledge of data bit, phase estimation is not feasible in this scheme. 

In terms of GPS applications, the problem of data modulation can be overcome by establishing a parallel 
(non-dynamic) link between the GPS satellites and a control ground receiver which also simultaneously 
receives the frequency translated version of the GPS receiver signals [13]. The data demodulation and 
estimation is then performed at the ground receiver. Once the data modulation is removed from the GPS 
receiver signal, then the problem reduces to simply estimating the phase, frequency, etc. of an unmodulated 
r-f carrier. This latter problem of course, has been extensively studied in the literature. See for example 
[9] for a very efficient method for simultaneous estimation of frequency and its derivative, and [12] for some 
similar techniques, including a third and fourth order Kalman filter, as applied to the GPS receiver problem. 

In this publication, we address the original (and more difficult) problem of estimating the signal pa- 
rameters from a data modulated sinusoidal carrier and propose a novel simultaneous estimation- detection 
scheme whose performance is very close to the estimation schemes proposed for the case of the unmodulated 
carrier. In fact, for the case of high dynamics under consideration, there is virtually no loss in terms of the 
required (P/No) due to data modulation. Apart from the fact that in GPS applications this does away with 
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the necessity of having a parallel direct satellite to ground receiver communication link, the proposed scheme 
is very important in many other similar situations including (modulated) signal acquisition from NASA deep 
space probes undergoing high dynamics, where such parallel links may not be feasible. 

In the literature [14, 15] there are several techniques of data detection. If the signal waveform is known 
precisely (and does not change from bit to bit) then data can be detected coherently by using matched 
filters or correlation receivers, irrespective of the actual waveform. If the waveform (carrier) has a constant 
known frequency, then either a coherent or differentially coherent detection may be employed depending 
upon whether or not the phase of the carrier is known. In a decision-directed version of these techniques, 
the carrier phase and/or frequency are estimated by a phase-locked loop technique and the data detector 
becomes part of the loop. It is clear that these techniques may not be feasible under low signal- to- noise 
ratio and the high dynamics considered here, in that the frequency may not be even nearly constant over the 
detection period and under such low loop SNR conditions it may not be possible to acquire the lock with 
data modulation present. 

In an alternative solution as in the Costas loop [15], the data detection problem is bypassed by a mul- 
tiplicative technique. However, such a loop also suffers in terms of loss of the SNR due to the multiplicative 
noise term which can be excessive for high loop filter bandwidths required and low received SNR. 

It is thus apparent that schemes which incorporate data detection in a loop which in turn is made 
dependent upon the acquisition and tracking of the loop may not be desirable under such high dynamics 
conditions since the loop may not acquire to start with (and may lose lock during tracking). Therefore, 
in this publication we propose and investigate techniques where data can be detected even under open 
loop condition. The fact that this information can subsequently be used to obtain an overall closed-loop 
configuration is only of secondary importance. 

The proposed scheme is based upon estimation of the received signal “local” (data dependent) param- 
eters over two consecutive bit periods, followed by a detection of a possible jump in these parameters. The 
presence of the detected jump signifies a data bit transition which is removed from the received signal. This 
effectively demodulated signal is then processed to provide the estimates of “global” (data independent) 
parameters of the signal which also provide a prediction for the “local” parameters over subsequent bit 
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periods. 

The estimation algorithm for both the “local” and “global” parameters estimation is an improved 
version of the Kalman filter. For “local” parameter estimation it is necessary to use an algorithm capable of 
estimating both the phase and frequency (data is phase modulated). However, for the “global” estimation 
algorithm it is not necessary to estimate the phase and thus the Kalman filter may be substituted by the 
frequency estimation algorithm of [9], resulting in a marginal reduction of the required P/N$. 

One of the key features of the proposed algorithm is the introduction of two different schemes which 
can provide an improvement of up to 3 dB over the Kalman filter (and phase or frequency locked loops) 
as applied to phase and frequency estimation under low signal- to-noise ratio (SNR) conditions. Both of 
these schemes are based upon the recognition of two available error signals with nearly independent noise. 
In conventional implementations, one of these error signals is simply ignored. Here we propose that by 
reprocessing the observations over an optimally selected period, one can exploit the other error signal as 
well. It may be emphasized that the improvement here is not due to better linearization as in iterated Kalman 
filter (where the improvement increases with the iteration interval) but rather due to the fact that if the 
reprocessing interval is optimal, the two sets of measurements are nearly independent, resulting in nearly 3 
dB improvement of performance under low SNR. In an alternative scheme, relevant to classical phase-locked 
loop structures, we propose an adaptive Hilbert transform technique resulting in similar improvements. 

In section 2 we present the signal model with the optimal solution, followed by the proposed simultaneous 
estimation and detection scheme in section 3. Section 4 presents the estimation algorithm in some details 
with the two proposed modifications. In section 5 we give simulation results where the algorithm is applied to 
the GPS receiver parameter estimation problem. The last section of the publication contains some concluding 
remarks. 
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2. THE SIGNAL MODEL 


Consider the problem of simultaneously estimating the high dynamic phase process 0(f) and detecting 
the data modulation from the measurements y(t),z(t) in (1) made in the presence of zero mean additive 
white Gaussian noise n,(f), n q (t). 


y(t) =ASin(0(f) + n D(t)) + t><(f) 


( 1 ) 


z(t) =ACos(0(f) -f vD(t)) + t^(f) ; 0 < t < To 

In (1) D(t) is a binary digital waveform taking possible values 0 or 1. Assume that over a sufficiently small 
estimation period, 


0(f) = #o + wot + To t 2 + ^ot 3 (2) 

for some unknown parameter vector tpo = [#o wo To ^o] ; . Thus {t/(f), ^(t)} represent the in-phase and 
quadrature components of a received signal s(t) and u 0 is the offset between the frequencies f c of the carrier 
reference signal r(f) and the received signal representing the Doppler, while the Doppler rate is equal to 2 to 
etc. In case of positioning systems, the phase derivative ©(f) is related to the physical velocity trajectory 
Vd(t) of the receiver as 


©(<) = 2 x^v d (t) (3) 

where c is the velocity of light. Similar normalization of the second and third derivatives of ©(f) yields the 
acceleration and jerk (second derivative of velocity) respectively of the physical trajectory. 

Sampling the inphase and quadrature signal waveforms y(t) and z(t) at a rate of l/T s one obtains 


y(k) = ASin(0(fc) + vD(k)) + t u(k) 


(4) 


z(k) = ACos(Q(k) + *D(k)) + »,(*) ; D(k) = 0, 1 

The signal y(t) (and similarly z(t)) may be depicted as in Figure 1. 

Optimum Solutions 

It is desired to obtain the estimate of parameter vector and detect data bits d lf . . . ,d n on the basis 

of measurements Y(M) t Z(M) “ {y(k),z(k);k = 1 where M=Nn, n is the number of received 
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digits and N is the number of samples/bit. The optimum solution can be obtained by maximizing the 

a-posteriori probability density function of V’O) rfi, • • • , d n conditioned upon Y(N), Z(N) and denoted by 

p(V>o,di, . . . , d n /Y(M), Z(M)), Such a maximization can be achieved by first computing the 2 n possible 

conditional probability density functions p, = max p(ip o|^(A0> Si) at their respective optimum 

*l>0 

point rpi. In the expression for p,-, Si denotes the ith possible binary sequence l , . . . , Si n for the received 
digits. This is followed by selecting the maximum out of these 2 n values. Thus if p m = max(pi, . . . ,p 2 n) for 
some m between 1 and 2 n , then rpQ Pt = and d \ y . . . , d n = S m . 

To evaluate the ith probability p f *, obtain the set of modified measurements Y(M),Z(M) defined as 


y{k) = Sijy(k) 

Y(k) = Sijz(k) ; (j - l)N + 1 < k < jN ; 
and then computing pi = max p(^o|^ (M ), Z(m)). 


( 5 ) 


When there is no a-priori knowledge about tpo, one may instead maximize p(Y(M), Z(M)\xp 0 ) with 
respect to ipQ. In case of Gaussian additive noise considered here, such a maximization is equivalent to 
minimization of the sum of prediction error squares, i.e., one simply minimizes 


M 

£{&(*) - -4Sin(0(Ar))] 2 + [z(k) - ^lCos(0(fc))] 2 } (6) 

k = 1 

with respect to rp 0 where Q(k) = ^[1 t * t\ with tk = kT 8 . The approximate minimization of (6) can be 
achieved by a recursive algorithm like the Kalman filter. 

In this publication we present a computationally simpler suboptimal version of the optimal estimator 
given above. Such a simplification is achieved by considering the case of n=l first and observing that one 
could estimate a modified parameter vector ip Q = tpo + Idiir with /' = [1 0 0 0], directly without an explicit 
knowledge of the data bit d\. The estimate of rp 0 in turn permits the computation of the total phase 
0(fc) = d\ 7r + Q(k) at the boundary of first bit (t — T&— ). A similar computation with the measurements in 
the second bit interval permits the computation of total phase d 2 7r-|-0(fc) at t = Ty f . Now it is obvious (see 
figure 1) that 0(fc) does not have any discontinuity at t = 7$ and thus any difference among the above two 
estimates can result either due to data transition or noise. Thus, if the signal-to-noise ratio is adequate, it 
is then possible to detect such a data transition with “small” probability of error. In this case the transition 
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can be eliminated from the measurements during the second bit period and the estimation of t/> 0 can be 
updated on the basis of the additional measurements. 

Assuming that the first bit is zero (without any restriction on the algorithm) the above decision-directed 
method effectively removes data modulation and provides at time k, the following modified measurement. 


y(k -N) = j4Sin(0(fc - TV) + nD(k - TV) + nD(k - TV)) ± v,(Jb - TV) 

( 7 ) 

z(k - TV) = ACos(Q(k - TV) + *D(k - TV) + itD{k - TV)) ± v q (k - TV) 
where N represents the delay equal to the number of samples in one bit period introduced in the physical 

realization and D(k) is an estimate of D(k). If D(Jc) = D(k) then of course, equation (7) reduces to 


y(k) = ASin(0o 4* wo T 8 k -f 7oT # 2 fc 2 4- fioTf k 3 ) H- 


( 8 ) 


z(k) = ACos(^o 4 u)oT s k 4 70 Tjk? 4 ^o^ a 3 ^ 3 ) 4 v q{k) 

From (8) one can easily estimate u> 0 , 70,^0 etc. for example, by an extended Kalman filter. The effect of any 
decision errors made in the detection of D(k) would be to introduce errors in the parameter estimates. But 
intuitively, such an effect would be most significant upon 6 0 , a parameter which is not significant in itself 
under the present application. 
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3. SIMULTANEOUS ESTIMATION AND DETECTION SCHEME 


In this section, we present a novel combined estimation/detection scheme which simultaneously detects 
data bits and obtains estimates of signal parameters such as carrier phase, frequency, etc. (related to receiver 
dynamics) in a sequential manner. The proposed scheme is recursive in both the number of data bits and 
the observations within any one data bit. The procedure described here effectively involves two mutually 
coupled estimation processes. In one of these estimation processes, we obtain the estimates of the signal 
parameters (phase, frequency, etc.) in the vicinity of possible data transitions (symbol boundaries) on the 
basis of measurements obtained within the current data bit. These estimates, which are dependent upon 
both the data and the receiver dynamics, and termed “local” parameter estimates, are then used to decide 
whether or not a data transition has actually occurred. On the basis of this information, data modulation 
is removed from the received signal and the modified observations are reprocessed to update the “global” 
parameters (dependent only upon the receiver dynamics and independent of data modulation) by taking into 
account the additional observations during the current “detected” bit. 

A block diagram of the proposed estimation/detection scheme is depicted in Figure 2. It may be 
noted that a closed-loop configuration involving the feedback correction signal to the reference numerically 
controlled oscillator (NCO) is required to keep the signal frequency 0(t) at the phase detector output within 
the bandwidth of the filter (accumulator) following the phase detector. Such a correction at a rate equal 
to a submultiple of bit rate consists of simply transferring the estimate of frequency (and possibly that 
of frequency derivative as well) to the NCO. In the following, we first describe in some detail the overall 
estimation/detection scheme without specific details of the estimation algorithm itself. In this paper we apply 
a modified Kalman filter but other appropriate recursive parameter estimation algorithms may be adapted 
to the proposed framework. The proposed estimation/detection scheme involves the following recursive (in 
number of data bits) steps. 

Recursion in Data Bits 

Step 1 : Let N be the total number of samples in any one bit period (assumed to be an integer) with the first 
sample appearing at time t = 0 + and the Nth sample occurring at time T$ — T s where 7$ denotes the bit 
period. With a recursive algorithm to be described subsequently, we obtain the estimates of the parameter 


7 


vector ip 0 = {So wo 7o ^o] T at time t = 0+ on the basis of first N measurements. Here #o = 0q 4- 7rdi with d\ 
denoting the first data bit with its value 1 or 0, i.e. denotes the total phase at t = 0 + . 

Step 2 : In order to detect bit transitions we need to obtain estimates of ‘local” parameters at the bit 
boundaries. Thus to detect the possible transition at the end of first data bit period, it is required to obtain 
estimates of the parameters at t = and t = T&-f . Noting that the parameters at t = T&— are related to 
the corresponding parameters at t = 0 + by the following linear transformation, 


m-) = 


[0(T b -)l 


rl Tj T fc 2 7f-i 


' *0 " 

«m-) 


0 1 2Tj 3T t 2 


Wo 

7(n _ ) 


0 0 1 3T 6 


7o 



.0 0 0 1 . 


UJ 


one could predict the estimate of parameter vector rp(T b ) via the same linear transformation, i.e.: 


( 9 ) 


= S(T t )j> 0 ; 

. 00 ) 

/rT , . mod,— 

rl> P (Tb—) — ^ ) + ^ 2 ") 

where S(Ti) denotes the matrix appearing in (9), /' = [1 0 0 0] and modulo 2 tt operation is performed on 
the first component of the argument. 

However, since we are interested in optimum modulo 2 it estimate of the phase ), a nonlinear 

function of ^ 0 , the modulo 2 tt version of ^(Xi— ) appearing in (10) and denoted a iJ>p(Ti,—) is not optimal 
even if ip 0 is optimum. Also in the subsequent kth bit period, the predicted estimate of the phase could 
differ by 7r due to data modulation if d* ^ do- Thus we simply set 0/>(T&— ) equal to its predicted value plus 
7r/2 or it may simply be set to zero. 

An optimum modulo 2ir estimate of tp(Tb—) is obtained by reprocessing the first N measurements but 
with their time and phase reference measured with respect to t = 7$ — as follows. With the initial estimate 
given by ^/>(Xi—), apply the recursive estimation algorithm to the set of following measurements to obtain 
an optimum filtered estimate of denoted tf F {Tb — ). 


y(j) = ASin(e L> i(j)) + t 

z(j) = ACos(0 L4 (j)) + v q (j) ; j = -AT, -~(N - 1) • • • , -1 


©l, i(j) = *(H-) + jT s u>(T>-) + 0X) 2 T (n-) + (;T,) 3 *(Ti-) (11) 
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We note that in (11), 0(7*— ) denotes modulo 2 tt phase at t — T b — and thus 0z,,i(j) may differ from 
Q(N + j) by an integer multiple of 2 tt. The fact that the error covariance of rp F (T b — ) may differ from the 
error covariance of rp P (T b — ) is apparent from the fact that the Cramer Rao error bounds are dependent on 
the time reference [1], 

Step 3 : assuming that all the parameters except 0 cannot change instantaneously, the predicted estimates 
of parameters at t = Ti-f (just after the possible data transition) are given by 


wp(7i+) = Q F (T b -) ; 7p(7*+) = 7 F (T b -) 


( 12 ) 


6p(T h +) = 6 F (T h -) ; 0 P (T b +) = 0 F (T b -) + w/2 

With the initial estimate given by tl> P (T b +) in (12), apply the recursive estimation algorithm to the set of 
following measurements {y(k), J(k) ; Jfc = 1, 2, • • • , N) to obtain an optimum filtered estimate of t/>(T b +) 
denoted xf> F (T b +). Under low SNR conditions 0 P (T b +) in (12) may be set to zero. 


y(j) = v(N+j) = ASin(Q Lt2 (j)) + Vi(j) 

z(j) - z(N + j) = ACos(0t >2 (i)) 4- ^(j) ; j = 0, 1 • • * , N - 1 


e £( 2 O') = Hn+) + }T,u(T b +) + (jT,) 2 7 (T 6 +) + (jT t ) 3 6(T b +) (13) 

In equations (11, 13) suffix L (Local) on 0 refers to the fact that the cumulative phase 0 is measured 
with reference to time T b — and T b + respectively. For computational simplicity, u> F (T b — ) etc. in (12) may 
be replaced by the corresponding components of ip P (T b —) of (10) and step 2, 3 may thus be performed 
concurrently. 

Step 4 : We note that irrespective of the dynamics involved (excluding the physically impossible case of 
instantaneous position change), Q(T b +) can differ from 0(T b — ) only if a data transition occurs at t = T b , i.e. 
d 2 ^ d\. In case of no transition we have that d 2 = di and F(7fc-f) = "G{T b — ). Thus we apply the following 
simple decision rule for data detection. 


d 2 = d r if \0(T h +) - 9(T b -)\ < tt/2 
/ d\ otherwise 

A detection error thus occurs only if the phase estimation error is greater than tt/2. 


(14) 
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Step 5 : In this step, the estimate of the global parameter vector V>o (parameters at t = 0 + ) is updated on the 
basis of modified measurements {i z\f(k) ; k = N, • * • , 2JV — 1}, where j/m (&) = 2/(k) if d 2 = di and 
yA/(ib) = — y(k) otherwise, with similar expression for z\f(k). Thus if there is no detection error made, one 
would obtain the same estimate for rpo as for the case of no data modulation. Any error in data detection 
simply appears as additional state noise in the algorithm. The significance of this step is that this yields 
(near optimum) estimation of the absolute phase 0 o (O + /2T6—), as against modulo phase, and from this 
absolute phase estimate 0 a bs(*/27& — ) (carrier phase in the absence of any data modulation) at any time t 
can be determined by simple linear transformation S(t) defined in (9). Equivalently using a dynamic model 
for phase [11], one could directly estimate total phase Q(k) at any time instance Ar. 

Whereas by necessity, 0(7},+) is based upon only the measurements {j/(fc), z(k) ; N < k < 2 N — 1} 
(previous measurements influence this only indirectly via w(T&-) etc.); the phase estimate 0 a bs (t/27&— ) is 
obtained on the basis of complete set of measurements up to instance 2 N — 1 and refers to the carrier phase 
in the absence of data modulation (assuming that the first bit is 0). This more accurate phase estimate can 
then be used to obtain a more accurate (coherent) detection of data if a delay on data detection but not on 
parameter estimation can be introduced. 

Step 6 : Steps 2-4 are now repeated for detection of subsequent data transitions and step 5 updates the 
estimate of rpo on the basis of additional measurements during consecutive bit intervals. Note that in the 
first step of the algorithm, one would have 

= S(kT„)t(0 + /kT b -) ; k = 1, 2 . . . 

The generalization of other steps of the algorithm is obvious. 


10 


4. ESTIMATION ALGORITHM (RECURSION IN SAMPLES) 

We now consider the problem of estimating the parameter vector xpo = [#o To <5o]' in the signal model 
(1-4) when the data modulation has been removed as described above. For this purpose, the signal model is 
written in the following state space form. 


x(k + 1) = $x(k) + w(k) 
Z(k) = h(x( ifc), k) + v(k) 


m = 



v(k) = 


Vi(k) 

Vq( k ). 


h(x(k),k) 


ASin(Q(k)) 

ACos(0(fc)) 


e(k) = j'(k)x(k) , •/'(*) = [! t k t\ t 3 k ] 


(15) 


Jj'(fc) = Po w 0 70 s 0 ] ; $ = I 

where I denotes the identity matrix and {v(fc)}, {id(A:)} are assumed to be zero mean white Gaussian noise 
processes with covariances given by 


E[lj( k)u'(k)] = Q , E[v(k)v'(k)] = R = 


r 0 
0 r 


(16) 


Model (15) is similar to those of [11, 16] except that here we treat the model state as the unknown param- 
eter vector rpo rather than the time-varying phase and frequency process. The two models are completely 
equivalent if there is no data modulation or in case of modulation, no detection errors are made. As shown 
in the appendix, the extended Kalman filter equations for the above model are given by 


x(k + l/k + 1) = $x{k/k) + K{k + 1) € (Jk + 1) (17) 

where the gain K(k+ 1) does not depend upon the parameter estimates. The updates for the gain K(k + 1) 
and the prediction error 6 (Ar -h 1) are given by equations (A 10, A 12) of the appendix. 

Improved Kalman Filter : 

In the following, we propose a novel (and simple) modification to the Kalman filter, which results in 
about 3 dB improvement in terms of signal-to-noise ratio. It is observed that the prediction error £ ( k ) in 
(17) may be written as 
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€ (fc) = .ASin(0(*)) + nj(k) 


where 


0(jfc) = Q(k) - Q(k) ; n/(ifc) = v,(fc)Cos(0(i)) - v 9 (k)Sin(0(k)) (18) 

Now, we define a related error term r](k) as 


»,(jfe) = y(Jfc)Sin(0(ib)) + z(fc)Cos(0(jb)) 

which may be rewritten as 

rj(k) = ACos(0(ib)) + nq(k) 


(19) 


ng(jfe) = Vi(k)SinQ(k) + v q (k)CosQ(k) 

From the observation that {vi(ib), v q (k) } u(k)} forms an independent sequence and the fact that 0(fc) is 
measurable with respect to the noise samples up to time k — 1 , it follows that 


£["/(*)] = 0 ; E[n Q (k)] = 0 , E[n)) = E[n% ] = E[n 2 Q ] = r , £[n,(*)n Q (*)] = 0 (20) 

Thus nj(k) and ng(jb) are independent. Strictly speaking {n/(Jb)} or {riQ(k)} is not an independent sequence 
as n/(ib) is correlated to nj(k — 1) via 0(Jb). However, such a correlation is negligible* and is generally ignored 
in Phase Locked Loop Theory [16]. 

Thus we observe that an independent observation is available and exploit it in the following manner to 
improve the estimate of x(k). 

Among several possible methods, we select an iterative method wherein the measurements z(k)} 

are divided into groups of M measurements with M equal to N or an integer submultiple of N. Each of 
these groups is then processed by the estimation algorithm 2 or more times. Considering the first such 
group with k between 1 and M, then under low signal-to-noise ratio (SNR) condition, the estimate of total 
* Essentially O(Ar) is negligible compared to Q(k) and thus Q(k) in the expressions for nj(fc), nq(k) may 
be replaced by ©(Ar) which is a process independent of the additive noise process (v^(Ar), 
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phase ©(&) in the first iteration denoted by ©*(£) would be considerably different from its estimate 0 2 (£) 
in the second iteration, i.e. |0 2 (Jb) — ©(Jb) | << |0 1 (ib) — ©(Jb)|. It is easily seen that the corresponding 
noise samples appearing in the prediction error term in (18) and denoted by n}(b),n 2 (b) respectively have 
a correlation coefficient equal to Cos(© 2 (fc) — © 1 (b)). Thus if the difference 0 2 (b) — ©*(&) is approximately 
equal to 7t/2, the corresponding noise samples are nearly independent and a signal processing gain of about 
3 dB is realized. This effectively corresponds to using a prediction error rj(k) given by (19), in the second 
recursion. For a given signal- to- noise ratio, the value of M can be optimized to achieve maximum possible 
improvement. Note however, that no such gain can be realized under high SNR case as under such conditions 
©^b) = ©(fc) — © 2 (b). These observations have been verified by computer simulations. It may be stated 
that improvement here is not due to iterative linearization of the nonlinearity as is usually the case but 
rather due to very specific properties of the noise n/(Jfc). 

Remarks : 1. As is usually the case, K(k+ 1) in (15) being independent of measurements may be precomputed 
and can possibly be replaced by its steady-state value. 

2. Here we have worked with a complex baseband signal. However, similar ideas can be applied to scalar 
sampling case [3, 11] where a Phase- Locked Loop structure results. 

Modified Phase Locked Loop with up to 3 dB Improvement 

In this section we also provide an explicit implementation in terms of a phase-locked loop (PLL), see [3, 
11] for details. Such an implementation is shown in Figure 3, where to keep consistency with the notation 
of previous sections, we represent both the received bandpass signal s(k + 1) and the NCO (numerically 
controlled oscillator) signal in terms of the same frequency lj c . Any difference can of course be absorbed in 
u>o, a part of the state vector. All other symbols in the figure are consistent with the rest of the publication 
except that ar(Ar) and $ would correspond to their use in [11]; i.e. x(k) here consists of instantaneous phase, 
frequency, etc. rather than the values of these at some fixed time instance. Note also that the bandpass 
noise v(k + 1) of Figure 3 has the following quadrature representation in terms of v f *(Jk) and v q (k) of (3). 

v(k) ~ Vi(k)Sin(u c t k ) + v q (k)Cos(u> c tk) 

The loop filter may be designed either on the basis of steady-state solution of the Riccati equation (Kalman 
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Gain) or via any other design procedure applicable to PLL filter design. 

To see the possible 3 dB improvement due to Hilbert transform, simply note that if 0((fc) is equal to 
kT s Q for some frequency error u> with rj(k) given by 


T](k) = ACos(u>kT t ) + ng(fc) 

then its Hilbert transform is simply 


( 21 ) 


fj(k) — ASin(QkT 9 ) + fiQ(k) 
= AS\n(Q(k)) + ng(Ar) 


( 22 ) 


Moreover, the independence of nj(k) and ng(Jb) also implies the independence of n/(fc) and nQ(k) due to 
Gaussian assumption on noise and the fact that the Hilbert transform is a deterministic linear transform. 
The fact that u) is actually time- varying is taken into account by making the Hilbert transform adaptive, i.e. 
in time-domain we replace the Hilbert transform by its following adaptive discrete time version; 




( 23 ) 


where A is some adaptation constant with 0 < A < 1. For A = 1, (23) corresponds to the discrete time 
version of the standard continuous-time convolution definition (24). 


w = (rb) v(r)dr (24) 

r 

The parameter A is selected so that over an interval T s /(l — A), the frequency u> remains nearly constant. 
Thus the selection of A depends upon the input signal dynamics, loop bandwidth, etc. It may be noted 
that the implementation of Hilbert transform by necessity requires a delay in the loop and thus the actual 
realized gain may be less than 3 dB depending upon the amount of such a delay. 
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5. SIMULATIONS 

In the following we present simulation results obtained when the algorithm is applied to the signal 
parameters acquisition and tracking for a high dynamic GPS receiver [7, 12]. For the purposes of this 
simulation we assume that the pseudo random code has been removed from the received signal and symbol 
timing has been acquired. Such a signal is then quadrature sampled at a rate that is an integer multiple N of 
the bit rate 1/T&, the latter assumed to be 50 bps. For the purposes of simulations, we consider two different 
sampling rates of 750 and 1000 samples/sec corresponding to the value of N equal to 15 and 20 samples/bit 
respectively. The simulation results are presented for different values of P/N 0l where P denotes the received 
carrier power and JVo is the one-sided power spectral density of the receiver bandpass noise. The parameters 
of equivalent baseband model (4) are then related as 


i 


A 2 1 PT b 
2 <r 2 N^N 0 ’ 


= E[vn 


(25) 


with -j^ (dB) = (dB) + 17. For the purposes of simulations we use normalized parameters as discussed 
in [8, 9]. The performance of the algorithm is completely invariant to such a simple transformation, but in 
implementation such a normalization can be very convenient as it keeps the condition number of the error 
covariance matrix small, and thus the algorithm would be more robust under finite precision arithmetic. 

The normalized sampling period T s n is selected to be .02 s with corresponding normalized parameters 
k>ojv, Jon and Son expressed in units of rad/sec. The parameters of interest i.e., the receiver velocity t ;<* in 
m/s and its derivatives in appropriate units can be expressed in terms of the above normalized parameters 
via 


v d = 


Nu>on 
2tt f e 

N^on c 


a = i)d = 2 
J = v d — 6 


2* Jc 

N 3 6on c 


(26) 


2 7T fc 

where c is the velocity of light and Sc is the received carrier nominal frequency equal to 1.575 GHz. Figures 
4-7 plot both the estimates of the normalized parameters as well as those of the receiver physical parameters 
t?d, a and J defined in (26). We consider first the acquisition problem where in the initial parameter vector 
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is some unknown constant. The modification to include the time- varying case is subsequently considered. 
For the simulation results presented below, a time-varying exponential data weighting factor given by 
(A7) is used with different values for the parameters Ao, A(0) and A max for the “global” parameter and 
“local” parameter estimation purposes. The values of these parameters used in simulations are specified 
below. 




Ao 

A(0) 

Amax 

‘Local” 

Estimator 

1.0 

.99 

.99 

‘Global’ 

* Estimator 

.995 

.99 

.998 


Note that for A 0 < 1, the value of A(Jfe) when computed according to (A7), it monotonically increases to 
1 as k — ► oo. In the simulations if A (Jfc 0 ) > A max for some k 0l we simply set A(Jfc) = A max for k > k 0 . 

Considering the estimation of “local” parameters, it is apparent that there is much more uncertainty 
in the phase parameter 0(7*- h) compared to other parameters (a good prediction of the other parameters is 
available from “global” parameter estimates). This fact was taken into account to marginally improve the 
performance by multiplying the (1, 1) element of the matrix \~ l (k)$P{k/k)& in (A5) by A(jb)/Xi with Xi 
= min (Ai, A(fc)). In the simulations Ai is equal to 0.9. 

A diagonal initial covariance matrix Pq{ 0/ — 1) (G denotes “global”) for the “global” estimation al- 
gorithm is selected with its diagonal elements equal to 10 6 , 10 6 , 10 3 , reflecting the initial uncertainty 
about the parameters. The corresponding initial covariances for the local estimations are simply given by 
S(mTi,)Po(mN/mN)S(mTi,) where S(m7*) is a transformation similar to 5(7*) of (6), and PoimN/mN) 
denotes the filter error covariance matrix for the “global” parameter estimator after mN measurements. A 
diagonal Q(k) matrix in (A5) is used with all its elements equal to zero except (1,1) element which is equal 
to 0.5 if k is an integer multiple of N in the global algorithm, or for k = 1 in the “local” algorithm to take 
into account phase noise due to data detection errors. Note that the extent of the “local” algorithms is 
confined to 1 bit of duration and thus the value of X(k) is reset to A(0) at the start of each new bit for fast 
convergence. 

Figures 4-7 present the results of typical simulation runs for the case when an unknown constant ac- 
celeration and unknown phase and frequency are introduced and a third order model for the state vector is 
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used. In the following we present results for u>j»(ib7i— ) and 7 jp(JbTj+) obtained in step 2 of the algorithm 
along with 0(kTb+) of step 3. Figures 4 and 5 present the simulation results for the case of 22.8 dB-Hz 
with the unknown acceleration of 50 g. Also, a randomly selected unknown phase of -2.3 rad. along with an 
unknown velocity error of about 25 m/s was included in the simulations. Such an error is expected when the 
algorithm is extended for tracking. In figure 4 we consider the case of no data detection (data bits assumed 
known). As is apparent from figures 4b and c, good estimates of both the acceleration and velocity are 
obtained within a period of about 10 data bits (200 ms). In this period the velocity error is reduced to about 
0.2 m/s with acceleration error of about 2 m/s 2 in various simulation runs. 

In figure 4c, the phase estimation error at the start of bit periods is plotted. Similar results are obtained 
for the phase estimates just before the data transitions at time ib7& — . From the drift in phase estimates as 
seen in the figure, it is clear that at the low CNR ratio considered here, phase coherence is maintained over 
short periods. However, such coherence is adequate for the differential detection discussed in the publication. 
In the simulation reported, six errors occurred over a period of about 40 bits. This may be compared with a 
DPSK detection (constant frequency) which results in a probability of error P e of about .06 at (PTb/No) of 
6 dB, and with coherent PSK detection (coherent phase) which results in P e of .01 [14]. Here, as the received 
signal frequency is time varying, the probability of error is even higher as may be expected. However, even 
under such a poor detection probability, good estimates of frequency and frequency rate can be achieved as 
shown in figure 5. 

In figure 5b, c, the estimates of velocity and acceleration are plotted for the case when data is detected 
simultaneously with parameter estimation. Compared to figure 4b, c the estimation errors are considerably 
higher. However, in a period of about 200 ms, velocity and acceleration estimation errors reduce to only 
about 1 m/s and 15 m/s 2 respectively, which are quite small. 

In figures 6 and 7 we present the estimation errors for the case when the sampling rate is equal to 20 
samples/bit. Keeping the sample SNR (A 7 /2a 2 ) the same as for figure 1, this corresponds to a (P/No) to 
24.1 dB-Hz. Comparison of figures 4 and 7 shows that the sampling rate does not significantly affect the 
estimation errors. Figure 7 presents the results for the case of 100 g step acceleration. 
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6. CONCLUSIONS 


We have presented a novel simultaneous estimation-detection scheme for the efficient estimation of 
phase, frequency and frequency derivative of a sinusoidal signal which may possibly be phase modulated by 
binary data and is received under noise. The proposed scheme is of a very general nature and can easily be 
adapted to situations other than those considered in the publication. For example, the scheme is equally 
applicable to situations where the data is QPSK on FSK (Frequency Shift Keying) modulated. The proposed 
algorithm is recursive in the number of data bits and within a bit period it is recursive in the number of 
signal samples. In terms of its application to GPS receivers, we observe from simulations, some of which 
have been reported here, that for a P/No of 23-24 dB-Hz, and with a 50-100 g acceleration, the scheme is 
capable of acquiring the velocity with an error of about 1 m/s in a 200 ms period. Thus, the algorithm 
can easily be adapted to high dynamic tracking situations as well. The most straightforward adaptation is 
simply to restart the algorithm every 200 ms. In fact, since the acceleration does not change by more than 
20 g (even for a 100 g/s jerk) over a 200 ms period, the actual performance could be expected to be better for 
tracking of 100 g/s jerk than reported here. Results of [7], which also consider data modulated carrier, show 
that an error of the order of 1 m/s is obtained for a P/N 0 of about 30 dB-Hz under similar dynamics. From 
this comparison it is apparent that the algorithm can achieve a very significant improvement of 6 dB over a 
previously proposed scheme based on Fast Fourier Transform (FFT) techniques. More extensive simulations 
of the proposed algorithm would provide a more precise measure of such an improvement. 
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APPENDIX 


Denoting by x(k/j) the estimate of x(Jfc) on the basis of observations up to time j, one obtains [11, 16] 
Filter State Equations: 


x(k + \/k + 1) = $x(k/k) + M(k + l)v(k + 1) 


(Al) 


v{k + 1) = z(k + 1) - h(x(k + 1/A), Jb + 1) 


(A 2) 


Filter Gain: 


M{k + 1) = P(k + l/k)h x (x(k + 1/A), k + l)*" 1 ^ + 1) (A3) 

Error Covariance Update: 

In (A3) P{k/j) denotes the error covariance matrix E[{x(k) — x(k/j)}{z(k) — x(k/j)} , ] y with a recursive 
update given by 


P(k + 1/A -I- 1) = P(A + 1/A) - P(A + 1/A )h' t <r-\k + 1)[P(A + l/k)h' x ]' (A4) 

P(A + 1 /k) = \~ 1 (k)$P(k/k)& + Q(k) (A5) 

c(k + 1) = h x P(k + \/k)h' x + R (A6) 

A(* + 1) = A 0 A(*) + (1-Ao) ; A(fc) < A max (A7) 


where A(£) is the exponential data weighting coefficient, Ao is some non-zero positive number less than 1, 
and h x denotes the derivative of the function h. 

To simplify the algorithm above for the given signal model, one obtains from (15), 
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h x (x(k/k 



ACos(Q{k)) ' 
-ASin(Q(k)) 


m 


e(i b) = j'(k)x(k/k - 1) 


Substitution of (A8) in (A6) results in 


a~\k) 


1 

r 2 *f ra(fc) 


a(ib)Sin 2 (0(A:)) -f r 
a(ib)Cos(0()b))Sin(©(ib)) 


a(ifc)Cos(0(ib))Sm(0(ib))' 
a(fc)Cos 2 (0(fc)) + r . 


a(jfe) = A 2 J'(k)P(k/k - l)J(k) 


Substitution of (A8) and (A9) in (A4) results in 


m 




P(k+l/k+l) = P(fc+l/jfc)-P(ib+l/Jfc)J(fc+l) J 4 2 [J'(fc+lW+lA)-/(fc+l)+'-]- 1 J'(k+l)P(k+l/k) (^110) 
The Kalman gain expression and the state update equations can similarly be simplified as 


M(k + 1) = yl(a(fc + 1) + r)~ l J(k + l)[Cos(©(fc + 1)) - Sin(0(* + 1))] (,411) 


x(k + l/k + 1) = <t>x(k/k) + K(k +!)€(*+!) 


K(k + 1) = [A 2 J'(k + 1)P(* + 1 /k)J(k + 1) + + l/k)AJ(k + 1) 

(- 412 ) 

G (k + 1) = y(k + l)Cos (Q(k + 1)) - z(k + l)Sin (Q(k + 1)) 

Equations (A12, (A5) and (A10) then constitute the extended Kalman filter equations for the signal model 


of interest. 
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Figure 2. Schematic Diagram of the Proposed Estimation/ 
Detection Scheme 
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Figure 4a. Errors in Estimates of Frequency at Bit Intervals 
(Acceleration 50 g) 
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Figure 4b. Errors in Estimates of Frequency Derivative at Bit Intervals 
(Acceleration 50 g) 
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Figure 5a. Errors in Estimates of Frequency at Bit Intervals with 
Data Detection (Acceleration 50 g) 
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Figure 5b. Errors in Estimates of Frequency Derivative at Bit Intervals 
with Data Detection (Acceleration 50 g) 
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Figure 6a. Errors in Estimates of Frequency at Bit Intervals 
(Acceleration 50 g) 
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Figure 6b. Errors in Estimates of Frequency Derivative 
(Acceleration 50 g) 


32 


ACCELERATION ESTIMATION ERROR (meter/sec2) 




ESTIMATION ERROR (rad/sec) 



Figure 7a. Errors in Estimates of Frequency at Bit Intervals 
(Acceleration 100 g) 


VELOCITY TRACKING ERROR (meter/sec) 



MAGNITUDE OF ESTIMATION ERROR FOR 
NORMALIZED PARAMETER Y (rad/sec 2 ) - 



Figure 7b. 


TIME IN NUMBER OF DATA BITS 


Errors in Estimates of Frequency Derivative at Bit Intervals 
(Acceleration 100 g) 


3 


ACCELERATION ESTIMATION ERROR (meter/sec2) 


TECHNICAL REPORT STANDARD TITLE PAGE 


1, Report No. 88-42 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Efficient Detection and Signal Parameter Estimation 
With Applications to High Dynamic GPS Receivers 

5. Report Date 

December 15, 1988 

6. Performing Organization Code 

7. Author (s) 
R. Kumar 

8. Performing Organization Report No, 

9. Performing Organization Name and Address 

JET PROPULSION LABORATORY 
California Institute of Technology 
4800 Oak Grove Drive 
Pasadena, California 91109 

10. Work Unit No. 

1 1 . Contract or Grant No . 

NAS7-918 

13. Type of Report and Period Covered 
JPL External Publication 

12. Sponsoring Agency Name and Address 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Washington, D.C. 20546 

14. Sponsoring Agency Code 
RE182 PX-644-11-00-0452 


15. Supplementary Notes 


. ract ijj^g publication presents a novel technique for simultaneously detecting 
data and estimating the parameters of a received carrier signal phase modulated by 
unknown data and experiencing very high Doppler, Doppler rate, etc. Such a 
situation arises, for example, in the case of Global Positioning Systems (GPS) where 
the signal parameters are directly related to the position, velocity and 
acceleration of the GPS receiver. 

The proposed scheme is based upon first estimating the received signal local 
(data dependent) parameters over two consecutive bit periods, followed by the 
detection of a possible jump in these parameters. The presence of a detected jump 
signifies a data transition which is then removed from the received signal. This 
effectively demodulated signal is then processed to provide the estimates of global 
(data independent) parameters of the signal related to the position, velocity, etc. 
of the receiver. 

One of the key features of the proposed algorithm is the introduction of two 
different schemes which can provide an improvement of up to 3 dB over ..the 
conventional implementation of Kalman filter as applied to phase and frequency 
estimation, under low to medium signal-to-noise ratio conditions. One scheme is 
based upon reprocessing (cycling) the measurements over an optimally selected 
interval while the alternative scheme proposes an adaptive Hilbert, transform 
technique. The overall complexity of the proposed algorithm is about three times 

the complexity of a single third order Kalman filter. 

17. Key Words (Selected by Author(s)) 18. Distribution Statement 

Spacecraft Communications, Command, and Unclassified — Unlimited 

Tracking; Communications; Navigation, 

Detection, and Countermeasures (General) j 
Navigation and Guidance 


19. Security Classif. (of this report) 

1 

20. Security Classif. (of this page) 

21 . No. of Pages 

22. Price 

Unclassified 

Uhclassified 

34 



JPL 0184 R 9183 






















